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Laboratory  Measurements  of  the  Water-Air  Flux  of  Dimethylsulfide  using  a  Wind  Wave  Tank 
J.  W.  H.  Dacey  and  N.  M.  Frew 

Estimating  the  flux  of  dimethylsulfide  (DMS)  from  the  ocean  to  the  atmosphere  is  an 
important  biogeochemical  problem  (Andreae  and  Crutzen  1997).  Sulfur  aerosols  are  an 
important  contributor  to  optical  haze  in  the  marine  atmosphere  as  well  as  a  possible  modulator  of 
global  climate.  Shaw  (1983,  1987)  first  suggested  that  atmospheric  sulfur  might  be  involved  in 
climate  regulation.  Charlson  et  al.  (1987)  proposed  that  volatile  organic  sulfur  emitted  from  the 
ocean  was  involved  in  the  formation  of  cloud  condensation  nuclei,  affecting  the  albedo  and 
radiation  balance  of  the  earth.  Charlson  et  al.  (1992)  estimated  that  the  direct  (backscattering  and 
reflection  of  solar  radiation  by  sulfate  aerosols)  and  indirect  (cloud  albedo)  effects  of  sulfate 
aerosols  might  counteract  the  climatic  forcing  by  trace  greenhouse  gases.  The  relative  roles  of 
greenhouse  gases  and  sulfate  aerosols  in  global  atmospheric  change  remain  a  topic  of  debate 
(Schwartz,  1988;  Caldeira,  1989;  Falkowski  etal.  1992;  Langner  etal.,  1992;  Kiehl  and 
Briegleb,  1993;  Bates  et  al.,  1997). 

Dimethylsulfide  (DMS)  is  the  most  abundant  volatile  sulfur  compound  in  seawater  and  is 
widely  accepted  as  dominating  the  flux  of  biogenic  sulfur  to  the  atmosphere:  it  constitutes  about 
half  of  the  global  biogenic  sulfur  flux  to  the  atmosphere  (Andreae,  1990;  Erickson  et  al.,  1990). 
The  concentrations  of  DMS  in  the  surface  mixed  layer  of  the  ocean  are  orders  of  magnitude 
higher  than  the  equilibrium  concentration  in  the  marine  atmospheric  boundary  layer.  Thus,  a 
strong  concentration  gradient  exists  to  drive  a  flux  of  DMS  into  the  atmosphere  (Andreae,  1986; 
Bates  et  al.,  1987).  The  sea-air  exchange  of  DMS  is  affected  by  a  number  of  physical  factors, 
including  seawater  temperature,  wind  speed,  surface  turbulence  and  gas  diffusivity  (Liss  and 
Merlivat,  1986;  Wanninkhof,  1992;  Saltzman  et  al.,  1993).  Saltzman  and  Cooper  (1989) 
concluded,  based  on  variations  in  DMS  concentrations  in  the  atmospheric  marine  boundary  layer, 
that  DMS  transfer  may  be  slower  than  predicted  by  the  standard  Liss-Merlivat  transfer  model.  A 
more  recent  assessment  by  Yvon  et  al.  (1996)  suggests  that  flux  may  be  higher  than  these  flux 
models  predict.  All  estimates  of  DMS  flux  have  depended  on  the  assumption  that  DMS  flux  is 
water-side  controlled,  although  the  relatively  high  solubility  of  DMS  suggests  that  its  behavior 
may  be  transitional  to  air-side  control. 

The  flux  of  gases  across  the  air-water  interface  is  understood  to  involve  three  transfer 
processes  (Liss  and  Slater,  1974):  (1)  transfer  from  the  bulk  aqueous  solution  to  the  surface  water 
film,  (2)  transfer  from  the  surface  water  film  to  the  surface  air  film,  and  (3)  transfer  from  the 
surface  air  film  to  the  bulk  atmosphere.  The  rate  of  gas  equilibration  at  the  interface  between  air 
and  water  films  is  taken  to  be  faster  than  transport  through  the  sublayers  on  either  side  of  the 
interface.  The  rate  of  transport  is  therefore  controlled  by  transport  to  or  from  the  interface 
through  the  viscous  sublayers  in  the  water  and  in  the  air.  The  relative  importance  of  the 
hydrodynamics  of  these  sublayers  is  determined  by  a  number  of  factors. 

The  rate  of  transfer  of  relatively  insoluble  permanent  gases  between  the  water  and  the 
atmosphere  is  controlled  by  the  flux  through  the  water-side  boundary  layer.  Mass  balance 
equations  indicate  that  higher  gas  solubility  increases  the  significance  of  the  air-side  control  of 
flux,  as  the  rate  of  equilibration  between  the  water  and  the  surface  air  film  has  a  disproportionate 
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impact  on  the  fugacity  on  the  air-side  boundary.  The  transfer  of  a  relatively  few  molecules 
quickly  diminishes  the  partial  pressure  gradient  across  the  interface.  The  fluxes  of  highly  soluble 
gases,  such  as  methylamines,  ammonia  and  water  vapor  itself  are  controlled  by  flux  through  he 
air-side  sublayer. 

Our  main  objective  in  this  project  were  to  parameterize  the  flux  kinetics  of  DMS  while 
simultaneously  measuring  water-air  transfer  of  02,  SF6,  Ne,  CH4,  and  He  under  varying  wind 
speeds  and  hydrodynamic  conditions.  The  experiments  were  designed  to  allow  investigation  of 
the  air-side  transfer  resistance  and  the  role  of  surface  films  in  regulating  DMS  flux. 

Experimental  estimates  of  k600  for  gases  measured  in  this  study  are  plotted  in  Figure  1. 

The  data  follow  the  non-linear  pattern  of  increasing  flux  with  increased  wind  speed  familiarized 
by  Liss  &  Merlivat  (1986).  Data  in  Figure  1  cannot  be  compared  directly  with  the  Liss-Merlivat 
plot  since  the  wind  speed  in  our  tank  experiment  is  measured  at  10  cm,  as  opposed  to  10  m 
elevation.  Our  data  are  revealing,  however,  in  illustrating  the  source  of  some  of  the  variability  in 
laboratory  and  field  measurements.  The  variability  in  flux  we  observed  under  similar  wind 
speeds  indicates  that  wind  speed  is  not  the  sole  factor  determining  flux  through  the  sublayers.  In 
the  case  of  water-side  controlled  gas  flux,  the  hydrodynamics  of  the  aqueous  sublayer  is  key,  and 
these  hydrodynamics  are  significantly  controlled  by  wind  stress,  ut.  In  the  field,  the  relationship 
between  u„  and  u  introduces  some  variability,  since  u*  depends  not  only  on  wind  speed,  but  also 
on  other  factors  such  as  sea  state.  This  does  not  explain  the  variability  in  the  wave  tank, 
however,  where  the  fetch  is  essentially  infinite  and  the  wave  field  is  established  quickly  and 
remains  constant  during  an  experiment. 
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Figure  1  Laboratory  wind-wave  tank  measurements  of  gas  transfer  velocity  (Sc=600) 
versus  wind  speed  (z=10  cm)  for  02,  SF6,  Ne,  CH4,  DMS,  and  He.  Data  represented  as 
open  symbols  and  solid  symbols  were  measured  in  fresh  water  and  seawater, 
respectively.  Fitted  lines  to  the  laboratory  measurements  of  Frew  et  al.  (1995)  are  shown 
for  comparison.  The  upper  bound  on  the  data  lies  near  the  clean  water  data  of  Frew  et  al. 
The  lower  bound  on  the  data  lies  near  experimental  data  with  surface  films  (0.03  uM 
Triton-X  100).  Standard  errors  are  plotted  for  each  point. 


Other  factors  control  the  relationship  between  w*  and  k600,  most  notably  variability  in  the 
physical  characteristics  of  the  water  surface  due  to  surfactants.  Small  changes  in  the 
concentration  of  surfactant  strongly  influence  the  hydrodynamics  of  the  aqueous  sublayer  (Frew 
et  al.,  1995).  This  variability  can  be  factored  out  through  the  dual-tracer  design  of  our 
measurements.  When  the  transfer  velocities  of  gases  are  normalized  against  the  transfer  velocity 
of  02  (all  k’s  corrected  for  diffusivity),  the  tight  functional  relationship  between  the  fluxes  of  the 
gases  becomes  apparent  (Figure  2).  This  relationship  indicates  that  much  of  the  variability  of 
flux  at  a  given  wind  speed  is  due  to  variability  in  the  hydrodynamic  forcing  of  flux.  Deviations 
from  the  linear  relationship  of  unity  slope  signify  other  processes  are  influencing  the  transfer  of  a 
gas.  Careful  examination  indicates  that  DMS  transfer  rates  lie  slightly  below  the  line  at  higher 
wind  speed,  meaning  that  DMS  behaves  somewhat  differently  than  the  permanent  gases  under 
the  experimental  temperatures  and  wind  speeds. 


Figure  2  Laboratory  wind-wave  tank  measurements  of  SFfi,  Ne,  CH4,  DMS,  and  He 
transfer  velocities  (Sc=600)  versus  the  transfer  velocity  of  02  (Sc=600).  Data  represented 
as  open  symbols  and  solid  symbols  were  measured  in  fresh  water  and  seawater, 
respectively.  Standard  errors  are  plotted  for  each  data  point.  The  solid  line  represents 
ideal  “water-side  controlled”  kinetics.  Dashed  lines  represent  predicted  km  for  DMS  at 
temperatures  0,  10  and  20C  using  our  formulation  of  ya. 

We  evaluated  the  importance  of  the  two  boundary  layers  by  considering  the  steady  state 
condition  in  which  the  mass  balance  for  gases  at  the  interface  is  given  by: 


when  the  flux  through  the  boundary  layers  of  the  air  phase  and  water  phase  are  equal  and 
respectively: 


ma"  =  Wa-CJV)  (2) 

<'  =  K(cs,w-cw) 


(3) 


where  ca  and  Cw  are  bulk  concentrations  in  the  air  and  water,  Csw  is  gas  concentration  at  the 
interface  between  the  water  surface  and  air  expressed  in  terms  of  aqueous  concentration,  and  p  is 
the  Bunsen  coefficient  (expressed  as  the  ratio  CJca  at  equilibrium).  The  transfer  velocities,  ka 
and  kw,  are  expressions  of  the  hydrodynamic  regime  in  the  air  and  water  boundary  layers.  We 
combine  eqns.  (1),  (2),  and  (3)  to  determine  the  relative  importance  of  the  air-surface 
concentration  difference  to  the  water-surface  concentration  difference  in  controlling  flux.  This 
ratio  becomes: 
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Using  the  fact  that  the  sum  of  the  air  fraction  and  water  fraction  is  unity  allows  equation  4 
to  be  arranged,  yielding  the  ratio  of  the  air  concentration  difference  to  the  total  concentration 
difference,  ya,  given  as: 


Pc<rc,,w  =  1 

K-cw  " 


(5) 


ya  is  the  air-side  gradient  fraction,  defined  as  the  gas  concentration  gradient  across  the  air- 
side  sublayer  layer  as  a  fraction  of  the  entire  gradient  from  the  bulk  water  to  the  bulk  air.  When 
is  small,  or  when  solubility  is  high,  ya  approaches  unity.  As  such,  therefore,  ya  is  a  measure  of 
the  error  in  assuming  the  flux  is  water-side  controlled.  To  date,  ka  and  kw  have  not  been 


U  (m/s) 

Figure  3  The  air-side  gradient  fraction,  ya,  for  DMS  as  a  function  of  wind  speed  at 
temperatures  0,  10,  and  20° C  in  seawater  are  shown.  ya  is  determined  from  p,  fcw,  and  fca 
as  shown  in  equation  5.  Values  of  p  come  from  Dacey  et  al,  1984.  Values  of  k w  and 
are  determined  from  tank  experiments.  The  value  of  is  determined  from  sparingly 
soluble  gases  and  the  value  of  is  determined  from  water  vapor  transport  (e.g.,  Kondo 
1975).  The  increase  in  ya  with  wind  speed  arises  because  &w  increases  faster  with 
increasing  wind  speed  than  ka.  Decreasing  temperature  raises  gas  solubility,  leading  to 
an  increase  in  ya 


quantified  simultaneously  in  oceanic  conditions.  In  situ  values  of  ka  (e.g.,  Kondo,  1975)  and  kw 
can  be  approximated  from  available  models. 

Due  to  the  moderate  solubility  of  DMS  (Dacey  et  al.,  1984),  ya  becomes  significant  at 
lower  water  temperatures.  Also,  increasing  wind  speed  has  a  greater  influence  on  the  aqueous 
sublayer  than  on  the  air-side  sublayer,  with  the  result  that  fca  becomes  increasingly  important  at 
higher  wind  speeds.  These  findings  are  summarized  in  Figure  3,  which  shows  the  influence  of 
temperature  and  wind  speed  on  ya.  We  conclude  that  air-side  control  may  be  a  significant  factor 
in  higher  latitudes  where  DMS  is  commonly  found  at  high  concentrations  in  surface  water 
(Berresheim,  1987),  and  where  both  low  temperatures  and  high  winds  are  typical.  As  ya 
approaches  0.4,  conventional  flux  models  based  on  wind  speed  parameterizations  may 
overestimate  DMS  flux  for  these  locales  by  60%. 
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